A simple second order differential 
equation, ovals and chaos 


In our youth as a consequence of our undying fascination with ovals we explored many means 
of generating them. In course of those explorations we experimentally arrived at a simple 
second order differential equation that generated oval patterns. It also taught us lessons on 
chaotic systems emerging from differential equations even before we actually explored the 
famous Lorenz and Rossler attractors for ourselves. The inspiration came from the very well- 
known harmonic oscillator which is one of the first differential equations you might study 
as a layperson: = —ax, where a is a positive constant that has some direct meaning 

in physics as the ratio of the restoration constant to the mass of the oscillator. With this 
foundation, we wondered what might happen if we used a cubic term instead which was in 
turn coupled to a periodic forcing from a regular harmonic oscillator. Thus, we arrived at 
the below equation with 3 parameters a, b, c: 
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It became obvious to us this equation has no solution in terms of elementary functions but 
we could semi-intuitively figure out that solutions are likely to generate oval-like figures. To 
test this out had actually solve it numerically. Just before this foray, a visiting relative, 
who was a college teacher, mentioned to us that the mark of man was to solve ordinary 
differential equations using the method of Ranga and Kutti on hand-held calculator. We 
had learned Euler’s method as a part of our education in basic numerical calculus but the 
this new method, apparently discovered by these south Indian savants, was something we 
had never heard of. Accordingly, we looked up a book featuring numerical analysis, which 
our father often recommended to us, and realized that the relative was actually referring to 
a generalization of the Eulerian method by the two sulapurusa-s Runge and Kutta. This 
differential equation offered us an opportunity to apply it to something interesting and we 
spent sometime writing a program to do the same. With that in hand, we could solve the 
above 2nd order DE by writing it as two linked 1st order ordinary DEs. 
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a = 6, b = 1.5 , c = 3.5 a = 4.5 , b = 3.5, c = 3.5 a = 6.5 , b = 3.5 , c = 3.5 a = 8.5 , b = 3.5 , c = 2 
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, c = 5.2 a = 4.75, b = 1.125, c = 5.2 a = 4.75 , b = 3.125 , c = 5.2 a = 1.75 , b = 8.125 , c = 5.2 
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Figure 1. Some solutions of §1 for different values of the parameters a, b, c with Xq = yo = 0. 

Consequently, we had an unfolding of the solutions for it which presented an interesting 
picture (Figure 1): We saw before our eyes the entire range of oscillations with simple pe¬ 
riodicity, quasi-periodicity, superimposed beats of different frequencies and various types of 
chaotic oscillations. Thus, we saw this simple DE recapitulate all manner of oscillatory phe¬ 
nomena we had encountered in nature: the changes is numbers of molecules in during the 
cell cycle, populations of organisms in the ecosystem, climatic patterns and the light-curves 
of variable stars. 
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If we plotted x, x we got our desired oval-like shapes. Notably, the cases where the solution 
is clearly chaotic we get a space-hlling entangling of near oval paths in this plot. 
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Figure 3. Plots of solutions with the b parameter changing while keeping a,c constant in each row. 


a = 3,b=2,c = 2 a=3,b=4,c=2 a=3,b=6,c=2 a=3,b=8,c=2 



Figure 4. x,x plots of the above solutions. 

From Figure 4 we can see that the amplitude parameter a and the wavelength parameter c of 
the sinusoidal term of §1 are the primary determinants of “shape” of the solution, while the 
parameter b associated with the cubic terms interacts with them to determine where chaos 
occurs. 
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